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D ■ Abstract 

Basing on the previously performed i^— matrix analysis of experimental data, 
we investigate, in the framework of the propagator matrix (D-matrix) technique, 
the 1100-1900 MeV mass region, where overlapping resonances /o(1300), /o(1500), 
<^ • /o fl530 ± 25o)) and /o(1780) are located. Neccessary elements of the D-matrix 

0^ . technique are developed. The Z)-matrix analysis confirms previous i^-matrix results: 

in the region 900-1900 MeV five scalar/isoscalar states are located. Four of them are 
fSj ' members of the two qq-nonets, while one state is an extra for the qq systematics, 

^D ■ being a good candidate for the lightest scalar glueball. The D-matrix analysis 

shows that this extra state, a candidate for the lightest scalar glueball, is dispersed, 
. ^ due to a mixing with qq- states, over three resonances: /o(1300), /o(1500), and 

p i' /o ( 1530 lb 250 ) • The broad resonance /o ( 1530 ± 259 ) is a descendant of the lightest 

p ' glueball carrying about 50% of the gluonium component, the rest of the gluonium 

5^ ! is shared between /o(1300) and /o(1500). 

> 
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1 Introduction 

At present time the mesons in the mass region 1000-2000 MeV are under intensive exper- 
imental studies. The problem these investigations face is the overlapping of resonances 
with identical quantum numbers: the widths of resonances in this region are compara- 
tively large, of the same order as mass differences of neighbouring resonances |l[]. This 
leads to an inevitable mixing of neighbouring resonances if they have common decay chan- 
nels. Consideration of this mixing is needed for a restoration of the initial-meson masses 
as well as for the determination of the quark content of the mesons. Definition of mixing 
is in particular important in the search for exotic mesons, such as glueballs and hybrids. 
Glueball (or hybrid) may be located in the vicinity of qq- mesons with identical quantum 
numbers that raises their mixing and leads to the dispersing of the glueball (or hybrid) 
component over several mesons. 

This paper is devoted to the calculation of scalar /isoscalar resonance mixing, IJ^'^ = 
00^^: experimental data prove that in this wave there is a state which is superflous for 
qq— systematics P, being possibly the lightest scalar glueball. The idea that in the mass 
region 1500-1700 MeV the lightest scalar glueball exists has been inspired by Lattice Glu- 
odynamics calculations |0, ^, ^. However, the quark degrees of freedom are not taken 
into account in these calculations properly, therefore they should be regarded as quahta- 
tive guide only: the gg-components can easily shift the glueball mass by 100-300 MeV. 
Glueball/gg-meson mixing was considered in refs. |^, 0, though without an implication of 
resonance decays into the mixing machinary. The K-matrix analysis demonstrates [0, § 
that the resonance decay processes may cause the mass shift by the same value 100-300 
MeV. The K-matrix amplitude, being unitary in the physical region and because of that 
treating the overlapping resonances correctly, does not reproduce left-hand side singular- 
ities of the amplitude which are related to the interaction of constituents. This means a 
neccessity to complement the K-matrix consideration by an analysis which restores the 
analyticity of the amplitude. 

Here we perform an analysis of the OO^^-wave in terms of the propagator matrix 
(or D-matrix). Corresponding technique, based on the dispertion relation N/D-method, 
reconstructs the amplitude with correct analytic properties. The detailed presentation of 
propagator matrix method is made in Chapter 2: we consider the mixing of overlapping 
resonances, the mass shifts caused by mixing and decay processes, and the decomposition 
of the final (i.e. physically observed) state in a series of initial states. Before only the 
00''"''-wave has been treated, however a generalization for other waves can be easily done, 
e.g. in the framework of the method given in ref. [^ |10 . 



Our D-matrix analysis of 00'*'^-wave is made in terms of gg-states, while K-matrix 
analysis of refs. p|, §] used the hadron language. When going from the hadronic K-matrix 
to the quark-antiquark propator matrix, the quark-hadron duality problem should be 
considered. It can be illustrated using standard quark potential model. In this model 
the gg-levels are defined by the potential which increases infinitely at large r, V{r) ~ ar 
(see Fig. 1). The infinitely increasing potential of confinement leads to the infinite set of 
stable gg levels: of course, this is a simplified picture of what is seen in the experiment. 
Only the lowest gg levels are stable with respect to hadronic decays. The higher gg 
states, 00"*"+ included, decay into hadronic channels: an excited (gg)a-state produces new 



qq pair, then the quarks {qq)a + QQ recombine into mesons which leave the confinement 
trap, forming a continuus meson spectrum at large distances. Conventional picture of the 
decay is shown in Fig. lb where the confinement interaction is depicted as certain barrier: 
the interaction at r < Rconfinement creates the gg-spectrum, while at r > Rconfinement a 
continuus spectrum of mesons takes place. In the experiment one can observe a continnus 
meson spectrum, so to extract information on qq spectra is a problem of 
(i) transition from hadron language to that of quarks and gluons; 
(ii) elimination of the infiuence of meson spectrum on gg- levels at r > Rconfinement- 

Sections 3 and 4 are devoted to the consideration of this problem. In Section 3 the 
necessary information is given concerning composite systems. In Section 4 the rules of 
quark combinatorics are summarized that allows to restore the quark content of a meson 
using the decays qq — meson — > two mesons. The same decays may serve as a signature of 
a glueball candidate. Section 5 is devoted to the calculation technique for the transition 
diagrams which are responsible for the mixing of levels created at r < Rconfinement- In 
Section 6 the discussed technique is applied to the analysis of 00"*"^ wave. Using the quark 
language for the transition diagrams, we investigate 1100-1900 MeV mass region where 
four OO^"*" states are located. Previously ||Tl|, the mass region 1100-1700 MeV was studied 
with resonances /o(1300), /o(1500), and /o(l530± ^^^j taken into consideration. The 
present investigation confirms the result of ref . 0] : the pure (gluodynamic) glueball is 
dispersed over three above-mentioned resonances, and the broad state /o f 1530 ± ^^A is 
a descendant of the glueball. The resonance /o(1780) has a small glueball admixture. 

2 D-Matrix Technique 

Here the D-matrix technique is presented in detail. First, we consider the propagator 
of non-stable particle (Section 2.1), then the propagator matrix for the two mixing reso- 
nances is constructed, followed by the generalization for an arbitrary number of resonances 
(Sections 2.2 and 2.3). Examples of a complete resonance overlapping are considered in 
Section 2.4: in this case one of the resonances created as a result of a mixing accumulates 
the widths of all other initial resonances. 

2.1 Propagator of Non-stable Particle 

Propagator of non-stable particle (resonace propagator) is given by the sum of diagrams 
shown in Fig. 2a, 2b, 2c, and so on. Scalar resonance propagator is equal to 

D{s) = {m^-s-B{s))-K (1) 

Here s = p"^, where p is the resonance four-momentum, and m is the input-state mass; 
d{s) = {m? — s)~^ is the propagator of an input state related to the diagram of Fig. 
2a, and B{s) is the loop diagram which describes the transition of input state into two 
particles. For the decay of scalar resonance into two scalar or pseudoscalar particles, the 



loop diagram B{s) may be written in a form of the dispersion integral as follows: 

Bis)= J '^A3V)m. . (2) 

J TT s' — s — iO 

fii and fi2 are meson masses in the loop diagram, g(s) is a vertex of the transition resonance 
-^ two mesons, and p{s) is invariant phase space of the two-meson state: 

p{s) = -^J[S - (/il + ^2^] [S - (/il - /i2)2] • (3) 

Complex resonance mass squared is determined by the propagator pole, m? — s — B{s) = 0. 
Below standard notation for it is used: s = M^ = M|. — zFM/j. Near s = M^, the real 
part of the loop diagram, B{s), can be expanded in a series over (s — M|.): 

D-\s) = {l + ReB'iM'^))iMl-s)-tg'is)pis) , (4) 



M^ 


= m^ - ReB{Ml) . 


Here we take into account that 




ImBis) = g'is)pis) , 


J 71 S' — S 



(5) 



The vertex g{s) has only left singularities of the scattering amplitude, thus being non- 
singular function in the resonance region. Smooth behaviour of g{s) justifies standard 
Breit-Wigner approximation for D{s), with resonance width being equal to: MnT = 
g\Ml)p{Ml)/{l + ReB\Ml)). 

For the resonance decay into several channels, one should make the following replace- 
ment in eq. (|ip: 

B{s)^Y.B''''\s). (6) 

n 

where index n refers to different channels. 

The loop diagram B^'^\s) is determined by eq. (|^), with a specification of vertices 
and masses: g{s) -^ fl'n(s) and /i, — *> pni- The scattering amplitude in the channel n is 
determined as 

A^is)=g^-\s)D{s)g^-\s). (7) 

The quantity {g^'^\M'^)y pn{Mf^) / Mr = Tn defines partial width of the resonance. 

2.2 Mixing of Two Resonances 

In the two-resonance case, the propagator of state 1 is determined by the diagrams of figs. 
2a, 2b, 2c type (transitions state 1^ state 1) and fig. 2d type (transitions state Instate 
2) . The sum of these diagrams gives: 

r, f ^ (2 u ( \ Bi2{s)B2i{s) y^ 



Here mi and m2 are masses of input states 1 and 2, while loop diagram Bab{s) is deter- 
mined by eq. (0), with the substitution g'^{s) — > ga{s)gh{s). Let the propagator matrix 
be introduced: 

Du Di2 

D21 D22 



D 



where the non-diagonal term Dab is the transiton propagator state b 
D-matrix is equal to: 



D 



1 



{M-^ - s){Mi - s) - B,2B2i 
The following notation is used here: 



Ml - s, 
B21, 



M( 



B12 

2 



(9) 
state a. The 



(10) 



M 



ra„ 



Baa{s) 



1,2. 



:ii^ 



Zeros of the denominator of propagator matrix provide the complex masses of resonances 
created as a result of the initial-state mixing: 

n(s) = {Ml - s){Ml -s)- BuB2i = . (12) 

Denote the complex masses of mixed states as M^ Mb- 

a) Model with constant Bgb 
As the first step, we consider a simple model: let us assume that loop diagrams Bab{s) 
depend weakly on s in the region s ~ M\ and s ~ M^. Let M^ and B12 be constants, 
then we have: 



M' 



A,B = liMl + Ml)±J^{Ml-Miy + BuB2i . 



(13) 



Small widths of the resonances 1 and 2 affect small imaginary part of Bi2. In this case 
eq. (^) gives the standard quantum mechanics result for splitting levels: a repulsion of 
mixed levels. 

The D-matrix can be written as a sum of pole terms, as follows: 



D 



Mf- 



+ 



M^s 



- cos sin a _i_ sin cos 6 



"Mj: 



"M|: 



- cos y sin a_ _i_ sin cos t 



Mi 



Mi 



sin^ 9 _|_ cos^ ( 



M 



Mi 



(14) 



where 



2n 1 1 

cos^ e = - + - 



\{Ml-Ml) 



2 2yi(M2-M|)2 + i?i2i?2i 
The states \A > and \B > are superpositions of initial states |1 > and |2 >: 

|y4 >= cos6'|l > -sin 6^12 > , 

15 >= sin6'|l > +cos6'|2 > 



(15) 



(16) 



b) General case In general case, when the s-dependence of loop diagrams is not neg- 
ligible and ImBij{s) is not small, the states \A > and \B > can be presented as a 
superposition of initial states as well. For this aim, let us consider the propagator matrix 



near s 



M' 



2. 



D 



Ii{s) 



Ml 



n'(Ml)(Ml-s) 



B2i{s) 

Mi (Ml) - 
B2i{M\ 



Buis) 

Mf{s)-s 



Ml 



Bi2{M\ 
MKMl) - 



Ml 



(17) 



Only singular (pole) terms are taken into account here. In the right-hand side of eq. f|l 
the matrix determinant is equal to zero due to eq. (O). Indeed: 



[M^iMl) - MI][MUmI) - Ml] - BuiMl)B2iiMl) = 
Eq. (12), n(M^) = 0, allows us to introduce the complex mixing angle: 



18) 



\A >= cos^aII > -sm9A\2 > 



(19) 



Then, the right-hand side of eq. ([T^) can be re- written, using 9a, as follows: 



D 



Na 



.Ml 



Ml 



COS^ 9 A 

sin 9a cos 9a 



— cos 9 a sin 9 a 
sin^ 9 a 



(20) 



where 



cos^ 9a 



Na 



Ml 



n'(Ml) 



[2Ml - Ml - Ml 



Ml 



2Mj - Ml - Mi 



2 ' 



sin^ 9, 



Ml 



Ml 



2Ml - Ml - Mi 



(21) 



Mj. In 



Remind that the functions Mf(s), Ml{s) and -812(5) in eq. (^) are fixed at s 
the case under consideration, when the angle 9a is a complex magnitude, the quantaties 
cos'^9a and sin^ ^^ do not relate to probabilities |1 > and |2 > in 1^4 >: the factors 
\/Na cos 9a and — ^/Na sin 9a are amplitudes for the transitions \A >^ |1 > and \A >^ 
|2 >. This means that corresponding probabilities are equal to | cos^^P and | sin6'Ap- 

To decompose the state \B > over initial states |1 > and |2 >, one should likewise use 
the propagator matrix near s = M^. Representing |i? > as 



\B >= sin6'B|l > +cos6'b|2 > , 
we have the following equation for D near s = M^: 

Nb 



D 



where 






Ml - s 



N. 



B 



n'(M|) 



sin^ 9b cos 9b sin 9b 

sin 9b cos 9b cos^ 9b 

2M| - Ml - Mf 



(22) 



(23) 



(24) 



cos^ 6b 



Ml - Ml 



sm^e 



B 



Ml - Ml 



2Ml - Ml -Mi' 2Ml - Ml - Ml ' 

In eq. (|2|) the functions Ml{s), Ml{s) and Bi2{s) are fixed at s = M|. 

When Bab{s) depends weakly on s and this s-dependence may be neglected, the an- 
gles 9a and 6b coinside. But in general case the angles are different. In this point the 
propagator matrix formulae differ essentially from those of standard quantum mechanics. 

Another characteristic feature of the discussed formulae is related to the mass shift 
of mixed levels: in standard quantum mechanics, we have a repulsion of levels associated 
with a conservation of the mean value, {Ei +E2)/'2 ( see also eq. ([TB|) ) . In the general case 
eq.([l^) can yield both the repulsion of mixed levels and their attraction. 

The scattering amplitude in the one-channel case is determined as 



Ms) = gais)Dabis)gbis) . 

For a multichannel case, one should introduce 

Babis)=j:Bi\s), 



(25) 



(26) 



where B^ is a loop diagram for the channel n with the phase space p„ and vertices g^\ 
gl"' . The partial scattering amplitude in the channel n is equal to 



A„{s) = g(-\s)D^b{s)gt\s) . 



(27) 



2.3 Propagator Matrix for an Arbitrary Number of Resonances 



Matrix elements Dab responsible for the transition of the initial state b into the state 
a, being a sum of diagrams of fig. 2 type, satisfy a set of linear equations: 



Dab = DacBcb{s){ml -s) 1 + 5ab{ 



m. 



:28i 



where Bab{s) is the loop diagram, and 5ab is Kroneker unit matrix. Let us introduce the 
diagonal propagator matrix d for initial (non-mixed) states: 



d 



{m\ - s) ^ 












m; 








[m% 



sY 



(29) 



Then eq.(p8D can be re- written in the matrix form, 

D = DBd + d 
therefor one has: 



D 



d'^~B 



(30) 
(31) 



The matrix d ^ is diagonal, so the inversed propagator matrix D ^ = {d ^ — B) is written 
in the following form: 



D 



-1 



-Bsiis) 



-Bu{s) 



^2 

-B. 



'32 (sj 



-Bisis) 

-B23{S) 

Ml-s 



(32) 



where Mf is determined by eq. (pT]). Inverting the matrix, one has for Dab{s): 



Dab{s) 



ba 



nW( 



(33) 



r(^-i)/ 



Here 11'^^^ (s) is the matrix determinant of D~^, and I^^ ^''(s) is the matrix supplement 
to the element [D~^]iia, i.e. the determinant of the matrix D'^ with the excluded 6th row 
and ath column. 

Let us write down n*^^''(s) as an example of the three-resonance case: 

U^'\s) = {M^ - s){M^ - s){Ml -s)- {Ml - s)B23Bs2 - {MJ - s)BsiBi3 - 

- {Ml - s)BuB2i - Bi2B2zBu - B1SB32B21 . (34) 

Zeros of n*-^''(s) give poles of the propagator matrix, which correspond to physical 
resonances created after the mixing of initial states. Denote the complex masses of reso- 
nances (or zeros of 11^^) (s)) as 



Ml, 



M|, 



M^,, 



(35) 



In the vicinity of the point s = M\, the pole terms provide the leading contribution. Ne- 
glecting the next-to- leading terms in eq. (|30|), one has a system of homogenious equations 
for Dab{s): 



D„ 



d-^ - B 



cb 







(36) 



The solution of this system, found with an accuracy up to common normalization factor, 
has a factorized form: 



f)iN) 



Na 



.Ml 



M2 



«?, 


aia2, 


aias 


02^1, 


"2 5 


a2a3 


"3"!, 


03^2, 


o^i 



(37) 



A^^ is a normalization factor which is chosen to satisfy: 



a\ + al + al + . . . + a]^ 



1. 



(3^ 



The complex coupling a^ is a normalized transition amplitude 



aa{resonance A — > state a) , 



(39) 



so the probability to find the state a in the resonance A is equal to: 

Wa = IttaP . (40) 

Analogous expansion of the propagator matrix can be done in the vicinity of other poles: 
d':,\s - Ml) = Nsj^^ , Di^\s ~ M^) = Ncj^^ . (41) 

The couplings satisfy the normalization conditions similar to eq.(38): 

f3f + (3l + ... + l3], = l, ^2^^2^... + ^^ = i. (42) 

But in general case, the completeness condition is not fulfilled for the inversed expansion: 

al + Pl + ^l + ...^l. (43) 

For the two- resonance case, this means that cos^ B^ + sin^ 0^ 7^ 1. Remind, however, 
that the right-hand side of eq.(|i3|) is equal to unity in the model with constant Bab'- see 
eqs. ([T3| ) - ([T6|) which correspond to the standard quntum mechnics formulae for mixed 
levels. 

2.4 Complete Overlapping of Resonances: Effect of Accumula- 
tion of Resonance Widths 

We consider here the examples which describe ideal situation with a mixing of completely 
overlapping resonances. These examples demonstrate in intact form the effect of width 
accumulation by one of resonances. 

) Complete overlapping of two resonances. 

Example 1: 

For the simplicity sake, consider the the weakly-s-dependent transition loop diagram, 
Bab] this case allows us to use eq. ([T3|). We suppose 



M2 = Ml- iMnT, , M2 = Ml- iMrT^ , (44) 

that means 

Eq. (|45|) can be valid even for positive gi and (72 if the integral over the region s' < M|. 
compensates the integral over s' > M|.. 
In this case 



i?i2(M^) = z^i(M^)(?2(M^,)p(M^) = tMn^T,T2 . (46) 

Substituting (H-dD into eq. ([ID, we have 



mXb = \{Ml + Ml) ± Ji(-zM2ri + iMlT,)^ + ( ^M^Vrir2 



Ml - zMr{T, + T2] 



(47) 



I Mr 

It is seen that, as a result of mixing, one state accumulates the widths of both initial 
states, Fyi = Fi + r2, while another one turns into a stable particle, Tb = 0. 

Example 2: 

Consider one more example when Re Mf and Re M2 are different but MiFi = M2F2, 
namely: 

M^ = Ml^ - iMT , M^ = Ml^ - iMV . (48) 



Then eq. ([T3|) reads: 



MIb = 2(^m + M'n2) - ^MT ± ^-(M|i - M^,)^ - M^T^. (49) 

This equation allows one to see the dynamics of poles with an increase of F. At 2MF <^ 
|M|^ — M|,2|, that corresponds to a suppressed mixing, one has two poles located near 
the positions given by eq. (0). With increasing F, the poles move to each other along 
the real axis. At 2MF = \Mj^^ — M^gl? the pole positions coinside with each other: 

MIb = ^{Ml, + Ml,) - tMT (50) 

With further increase of F, the poles move along the imaginary axis: as a result, there 
are two poles, one above another. At \Ml^ — M|,2| -C 2MF, one state is almost stable 
while the width of another resonance is close to 2F. 

b) Complete overlapping of three resonances 

Consider the equation 

n(=^)(s) = o (51) 

in the same approximation as in Example 1. Remind, n*^^-' is determined by eq. (|3^) . 
Correspondingly, we put 

Ml = Ml-s- iMnVa = x-i'ja, (52) 

ReBabiM],) = , (a^b). 
Here x = M^-s and MuTa = 7a. Then, with BatBba = i'^'yalb and B12B23B31 = 2^717273, 



eq. (0) reads: 

x^ + x^{i-fi + i-f2 + i-f3) =0. (53) 

So, the propagator poles are at 

Ml = Ml - zMh(Fi + F2 + F3), (54) 

Ml = M^ = Ml. 

Resonance A accumulates the widths of all initital resonances, while the states B and C 
turn into being stable and degenerate. Remark, the degenaracy goes off, when a weak 
s— dependence of loop diagrams, Bab, is taken into account. 
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3 Composite systems 

Up to now, the resonance is treated as a mixture of input states, and the loop diagrams 
Bab play crucial role in this mixing. In addition, they take another function, for they 
determine the compositeness of a particle. 

3.1 Entirely-composite particle 

Remind, following refs.f^ and [|12|, the main features of the N/D-method in a description 
of composite particles, like deuteron: this gives us an example of entirely-composite sys- 
tem. Entirely-composite particle appeares as a result of the interaction of constituents 
and reveals itself as a pole of the scattering amphtude of constituents. Correspnding 
scattering amplitude for multichannel case is 



glis) 



l-Bis) 



(55) 



where B{s) is a sum of loop diagrams given by eq. (P). This amplitude differs from that 
considered in Section 2 by the replacement {tti^ — s) —>■ 1: the absence of the input-particle 
propagator in the scattering amplitude is a signature of an entirely composite-state. The 
mass of the composite state, M, is determined by the condition: 

1 - B{M^) = 0. (56) 

Suppose that the mass M is located below all the thresholds, so we deal with a stable 
composite particle. In the vicinity of s = M^, the scattering amplitude can be represented 
as a pole diagram of fig. 3d: 



l-B{s) B'{^P){M^-s) M^-s' 
The vertex Gn{s) describes the transition composite particle —>■ n — channel constitents 

Gn{s) = ^MfL. (58 



The wave function of an entirely-composite particle is the N-dimensional Fock column as 
follows: 

Ms) 



^7v(s; 

where the partial wave function is determined as 



(59) 



Ms) = ^ . (60) 
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Definition of tlie partial wave function proceeds from the calculation of the composite- 
particle electromagnetic form factor, given by the sum of triange diagrams (fig. 3e), 
different channels taken into account in the intermediate state (for detail, see refs. P, p!2|): 

A = s^ + s'2 + g^ - 2ss' - 2sq^ - 2sV • 

Here, to avoid cumbersome expressions, we put the masses of constituents in the inter- 
mediate state equal to each other, /ii^ = /i2n = /in- The function A„(s, s', q^) in the limit 
g^ ^ satisfies the equation: 

A„(s, s\ g' -^ 0) = T^pn{s)5{s - s') , (62) 

This leads to the following relation: 

m-E/^^M.)^. (.3) 



no) ^E/v^-W 1^1^-1- («) 



Substituting ( pSj ) into right-hand side of eq. (p3|) , one sees that -F(O) is equal to unity: 

ds , , , Gnis) 

In terms of the relative momenta of constituents, eq. ([6^ ) reads: 
where k\ = 5/4 — /i^. The quantity 

^ ^n(^) , (66) 

is the probability for the composite particle to be in the state n, so 

l = EW^n- (67) 

n 

This means that the composite particle is entirely built of the constituents. 
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3.2 Non-Entirely Composite Particle 

Now return to the case when the amphtude is determined by a set of diagrams of fig. 
2 but with a pole of the amphtude located below all the thresholds. This means that 
we deal with a stable particle. In this case the Fock column which describes the wave 
function of the composite particle contains additional component, ipo, which corresponds 
to the point-like input state with bare mass m: 



^ 



^0 

Ms) 
Ms) 

iPn{s) 



(68) 



Consideration oiipn{s) is made in a way similar to previous case. The scattering amplitude 
for multichannel case is equal to 



glis) 



IV? — s — B(s) 



(69) 



Remind that B{s) is given by eq. @. The mass of the composite state, M, is determined 
by the following condition: 

Let us stress once again that M is supposed to be below all the thresholds, so this is a 
stable composite particle. In the vicinity of s = M^, the scattering amphtude can be 
represented as 



9lis) 



9\s) 



Glis] 



m?-s-B{s) (l + 5'(M2))(M2-s) M'^ - s' 



(71) 



The vertex (^^(s), responsible for the transition composite particle—* n-channel con 
stituents is equal to: 

gn{s) 



Gn{s) = - 

^l + B'{M^) 

Partial wave function in the channel n reads: 

GJs) 



i'ni.s) 



s-M^ 



(72) 



(73) 



As before, the definition of the partial wave function proceeds from the calculation of 
electromagnetic form factor: 



N "^ 



ds ds' G„(s) A„(s, s', g^)G'„(s' 



TT" 



{s - M^){s' - M^) 



(74) 



A„(s, s',g^) is determined by eq. (|6T|). Here, as compared to eq. (|6TD , an additional 
constant term should be included into the form factor, Wq, which corresponds to the 
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interaction of photon with a point-hke component (input particle with mass m): Wq 
is the probabihty for a particle to be in the input state. This leads to the following 
normalization condition: 



N °° 









W„.M^-f^.€i-s). ill) 



° 1 + 5'(M2) 
So, one has: 

The probability for the considered particle to be in the state n is equal to 

S;(M2) _ r dH 
1 + B'{M^) ~ J (2^)3 

N 

The magnitude J2 ^n is the probability for a particle to be in the composite mode. 

n=l 

3.3 Non-Stable Particles: Small and Large Distances in Transi- 
tion Diagrams 

Above we have discussed the stable composite particle. In the language of a potential, 
this corresponds to a level inside potential well: see fig. 4a. Now consider the case of 
non-stable particle (resonance): corresponding shape of the potential well and the level 
location is illustrated by fig. 4b. The level is mostly created by the potential at r < rg 
(where tq corresponds to the potential barrier maximum), while at r ^ tq the wave 
function describes outgoing particles, which are the resonance decay products. One may 
introduce the notions "inside the well" and "outside the well" in the simplest way, using 
step-function 0(ro — t), then the probaility for constituents to be bound is determined by 

Wn = j dh-4>l{r)Q{r^-r) . (78) 

Actually, we have introduced here the wave function 

CW=0nW0(^o-r), (79) 



which is equal to ipni''^) at r < ro and to zero at r > rg. Likewise, we can introduce the 
wave function for constituents at r > tq: 

0r*W = 0n(r)e(r-ro). (80) 

In the momentum representation, the wave functions read: 

^Wo«t(s) = f d^r(l)^/°^\r)e'^-'' , (81) 
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where k = |(/ci — /C2) and f = fi — r2. 

Presentation of the composite-particle wave function as a sum of short- and long-range 
pieces is the basic point in the analysis of O"*""*" wave. The gg-resonance being a composite 
state of quarks is determined by the diagrams of the type shown in fig. 5a: quarks interact 
by gluon exchanges. The quark interaction at r < Rconfinement results in the creation of 
bound states: in our analysis we approximate the small-r amplitude as a set of pole terms 
(see fig. 5b): 

E -^ ■ (82) 

However, the sum of pole terms does not take into account the interaction outside the 
confinement barrier, at r > Rconfinement- The quark loop with r > Rconfinement (fig- 
5c) should be saturated by a set of hadron loop diagrams (figs. 5d, 5e, and so on), in 
consequence with quark-hadron duality. Therefore, the inclusion of hadron loop diagrams 
into consideration is a neccessary step in accounting for the large-r quark interaction. 

In terms of the ^-representation, the separation of small-r and large-r contributions 
occurs due to a choice of vertices ga{s): suppression of the large-s region means suppression 
of the small-r domain, and vice versa. Of course, the use of the step function Q{r — rQ) for 
a separation of the large-r contribution is a simplification, which leads to the violation of 
analytic properties of the transition loop diagrams, B{s), in the physical region: actually, 
the small-r region switch-off should be done in a softer way. To keep the analyticity of 
i3(s)'s, one should care for analytic properties of vertices, ga{s). 

4 The Couplings: Rules of 1/A^- Expansion and Quark 
Combinatoric Relations 

When analysing experimental data, we approximate the amplitude at r > Rconfinement 
by the set of pole terms, see eq. (|5^: these poles are determined by the quark-gluon 
interaction. Quark-gluon content of the states related to poles reveals itself in coupling 
constants which govern the decay of these states into mesons, thus allowing to find out 
from the experiment the quark-gluon composition of input states. 

In this Section, first, the order of coupling constant magnitude is evaluated within 
the rule of 1/A^ expansion |T^ (Section 4.1). Then, the decay constants are calculated 
(Section 4.2) using quark combinatoric rules. These latter were used rather long ago for 
the calculation of secondaries in hadron-hadron collisions at high energies WM and in the 



J/\E' decay |T^. Calculation of the decay glueball -^ two pseudoscalar mesons within 



the quark combinatorics was carried out in refs. |]^, R, |TB[. In the same guideline, the 
decay couplings are calculated for the process qq — meson -^ two pseudoscalar mesons 
in Section 4.3. 
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4.1 Estimation of the Decay Couplings in the Framework of 
1/A^-Expansion Rules 

The rules of the 1/iV-expansion [|l^] (A^ = Nc = Nf where Nc and Nf are numbers of 
colours and light flavours) can be used as a guide in the soft quark-gluon physics. Let 
us estimate, in terms of 1/A^-expansion, the order of magnitude of the decay coupling 
constants Gluehall —>■ two mesons and Glueball —>■ qq denoted below as gG^mm and 

9 G^qq- 

For this aim, consider the gluon loop diagram in the glueball propagator, fig. 6a, 
which is determined by the coupling Gluehall — > two gluons denoted as Qc^gg- This loop 
diagram is of the order of unity if the glueball is a two-gluon composite system. From the 
other side, this loop diagram is proportional to QQ^ggN^ {N^ is a number of colour states 
in a loop at large Nc) , therefore 

B{G ^gg^G)r^ 9l^ggK ~ 1 , (83) 

or gc-^gg ~ 1/A^c • 
The coupling g c^qq is determined by the diagram of fig. 6b. Similar estimation leads to 

g G^gg ~ g G^ggglcoNc ~ 1/A^c • (84) 



Here gqcn is gluon-quark coupling which is of the order of Jl/Nc, while the factor Nc is 
related to the number of colour states in the triangle diagram. Then the process shown 
in fig. 6c determines the coupling of the Glueball — * two mesons decay in the leading 
terms of the 1/iVc-expansion: 



*mm 



gG^qqgl^gqNc-^l/Nc. (85) 



Here we take into account that the coupling gm^qq is of the order of l/\/Wc\ it is because 
the quark loop diagram in the meson propagator B{m ^ qq ^ m) (see fig.6d) is of the 
order of unity: 

B{m ^ qq^m) ^ 9m^qqNc ~ 1 ■ (86) 

The diagram of fig. 6e provides the contribution which is of the next-to-leading order 
of magnitude: 

altmm ~ 9 G-.gg9'QCD9l^qqN! ~ l/N^ (87) 

Low-lying glueballs may decay into the following two-pseudoscalar channels: 

n = TTTT, KK, rji], rji]', rj'r]' . (88) 

The ratios of couplings for the processes of fig. 6c can be calculated in the framework of 
the quark combinatoric rules. The important point here is the flavour symmetry violation 
in the production of the light quark pairs: uu, dd and ss. The production probabilities 
refer as 

uu : dd : ss = 1 : 1 : \, (89) 
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with A = 0.45 ± 0.05 [T^, [T8|. The couphngs for the glueball decay into two pseudoscalar 
mesons calculated within quark combinatorics are presented in Table 1 for both leading 
and next-to-leading terms, Qq^^^ and Qq^^^. Unknown dynamics of the decay is hidden 
into parameters G^ and G^^ . The decay coupling to the channel n is equal to the sum 
of both contributions: 

9G^mm{n) + gctmrnin) (90) 

The second term is suppressed, compared to the first one, by the factor l/Nc'. an experi- 
ence with calculation of quark diagrams gives us evidence that it leads to a suppression 
by a factor about 1/10. The sums of coupling constants squared obey the sum rules: 

Y.i9h-.mm{n)ri{n) = lGl{2 + Xr, (91) 

n ^ 

T.i9Gt^m{n)r I{n) = \gI,{2 + A)^ 

where I{n) is an identity factor for the particles produced (see Table 1). These sum rules 
reflect the fact that the decays under consideration are resulted from the cutting of the 
two-loop diagrams shown in figs. 6f and 6g, respectively: each quark loop contains the 
factor (2 + A), see eq. (|S|). 

4.2 Quark Combinatoric Calculations of the gg-Meson Decay 
Couplings 

Quark combinatoric rules can be applied to the (gg)a-meson decay if we know the quark 
content of the meson. Let the (gg)a-meson be the following mixture of non-strange and 
strange quarks: 

{qq)a = nn cos^ + nn sin^ , (92) 

where nn = [uu + dd)/\/2. Decay couplings of this meson to two pseudoscalar meson 
channels, eq. (|88|) , are determined by the diagrams of fig. 7a and 7b types. 

The process of fig. 7a provides the leading contribution, in terms of l/A'^c-expansion: 



ym{a)—*'tnm ym^qq c fTj" V / 

Similarly, one has for the process of fig. 7b: 

Jm[a)—tmm ^m^qq c ^m—tqq tit fTj" \ / 



The couplings of the decays m{a) —>■ tttt, KK, rjr], rji]', rj'r]' are presented in Table 2 both 
for leading and next-to-leading terms: g^ and g'^^ are parameters in which the unknown 
dynamics of the decay process is hidden. The decay coupling to channel tt, is a sum of the 
leading and next-to-leading terms, 

gm{a)^-mm\P') ~^ gm{a)—immv'') • \^'^) 
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These two terms are coupling constants in the most general form, thus giving all the 
varieties of decays. Examples with a special choice of coupling constants may be found in 
refs. [16,18], where it was suggested that g^k)^mm i^ ^^^ small as compared to gln(a)^mm- 
Let us emphasize a very important feature of the meson decay couplings, gm{a)^mm'- 
at $ = ^Giuebaii where 

tan^Glueball = Y 2 ' (^^) 

the ratios of meson couplings coincide with those for the glueball both for leading and 
next-to-leading terms. Namely, 

[9m{a)~*mm\''^)\'i>='!>Glueball ~^ 9G~*mmv'') ' l^' J 

It is due to the two-stage sructure of the glueball decay: on the first stage, the gg-pair 
is produced with the flavour content given by eq. (^91) that corresponds to $ given by 
eq. (pGl). The final stage is a decay of this gg-state. An important consequence of the 
two-stage structure of the glueball decay reads: basing on the ratio of the decay couplings 
to hadron channels only it is impossible to find out either we deal with the glueball or 



gg-state which has nn/ss content close to that of eq. (|89|) 



5 Resonance Mixing in Terms of the Quark Transi- 
tion Diagrams 

Quark-hadron duality means that one can analyse the OO^^-wave using quark or hadron 
languages. For our analysis this means the choice of a language for the description of loop 
diagrams — the transitions may be described using either quark or hadron states. In this 
Section the presentation of loop diagrams, Bab{s), is made in terms of quark states. For 
the description of quark states, the light cone variables are often used, so, as the first step, 
in Section 5.1 the loop diagram is written in these variables. As was mentioned above, the 
quark model deals with wave functions in the r-representation; in Section 5.2, Bab{s) is 
written using not vertices but wave functions. Simple wave functions with one-parameter 
dependence are presented in Section 5.3. 

5.1 Loop Diagram as a Function of Light Cone Variables 

Consider the loop diagram of eq. (^ using light cone variables. For this aim we write the 
phase space factor p(s) in general form: 

pW ^ mPMM - -^i2.mP-k. - fe)^5-^^5-^ (98) 

Rewriting B{s) with light cone variables can be done imposing Pz —>■ oo. Then, denoting 

kiz/Pz = Xi , rn? + k\ = m]_, (99) 
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one has for the right-hand side of eq. (|98D : 
1 dxidx2 



4(27r)2 X1X2 



(5(1 - xi - X2)d ki±d k2±6{P± - ki± - k2±) 



mj_L m2_L , 



<5(s + Pi-^^^-^^^). (100) 

xi 0:2 



Putting Pi = 0, we rewrite B{s) in terms of the hght cone variables: 



^(^) = 77^ / ^^^6{1 - xi - X2)d'k,^d'k2^6{h± + 4±)^^-^ . (101) 

4(^^7rj J XxX2 "S — s — zU 

Here s' = mlj_/xi + m2^/x2- 

5.2 Loop Diagram in the r-Representation 

To present the loop diagram as an integral in the coordinate space, let us rewrite it in 

the cm. ^-representation, where k = {ki — k2)/2 and ki + k2 = 0. Then, the equation for 

Bab{s) reads: 

f d^k fcio + k2o ga{s')gb{s') 
Bab[s) = / 7-—^ — — — , (102) 

s' = {kio + k2o)^ , kio = \Jmf + P . 
In terms of wave functions of the states a and b, it is equal to: 

o / N f d^k , , ,, {s' -ml){s' -ml) , , ,, ._^^. 

Babis) = J — -^ Ms') ., ,r, Ms') ■ (103) 

Here 



kio + k2o gai-s) 

In the coordinate representation one has for Bab{s): 

Babis) = fd\d'r'Mr) v^'^'Hlf- r^|; s)0,(r') , (105) 



^H)n.1. s) - f -^ e^-^' {s'-ml){s'-ml) 
^' '' ^~i (27r)3'^ s'-s-zO 

The wave function in the coordinate representation is determined as usually: 

d^k 



W 



e'^-'Ms) . (106) 



Constituent quark model makes it possible to estimate quantitatively the ^P^qq wave 
functions: making use of them, eqs. ( |103| ) and (|105|) allow us to calculate the transition 
loop diagrams. 
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5.3 Quark Transition Diagrams 

As was stressed above, when the resonance mixing has been analysed, two alternative 
languages can be used: the hadron language or the quark one. An advantage of the quark 
language is that the results of the quark phenomenology may be used as a quantitative 
guide. 

The quark loop diagram for the transition scalar state b -^ scalar state a is equal to: 



oo 



B„,(,, = / ^ ^-'f)"''')"''"'' TrlCk + m)(P - k - m)] = (107) 

J TC S — S — l[) 

TC s' — s — iO 

Am? 

where m is quark mass. In the trace calculations, we put P^ = s' and k"^ = m^: there 
is no momentum conservation in the dispersion representation of the loop diagram, but 
intermediate particles are on mass shell (for more detail, see refs. |1^, |19|). 

In Section 3, a receipt is given for the construction of composite-particle wave function. 
As a basis, one can use -B^^(M^) where Ma is the composite-particle mass: this quantity 
is equal, up to normalization factor (, to the wave function squared integrated over d^k. 
Applying this receipt to eq. ( |107| ), we get: 



BUM!) = / g;l ,.4f-ig) l' (S" ■ *)(S ■ *) - (108) 



</(&*°'*>*"'*> 



Here C = B'^(M2) and 



GJs) - »"'■') 



In eq.( |109| ) we introduce the spin-1 operator S. In standard representation it reads: 

5o= (0,0,1) 5±i = (i=,±-^,0) (110) 

The wave function 

■^a{k) = iS-k)Mk) (111) 

can be rewritten in the r-representation as follows: 

$„(r) = {S-k)Ur) = / (0 ^"^''^'^(^) • (112) 
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Here 

(S■r)Mr)=^{S■-^)Va{r), (113) 

Examples 

Now we present simple examples of the N^ Poqq-wave functions ( A^ = 1, 2), with a one 
parameter dependence. We suppose an exponential decrease of wave functions at large r. 
Then, the wave function of the A^ = 1 state in the r-representation reads: 

o n 

Parameter a is related to the mean radius squared < r^ > as follows: 

< r^ >i^Po= ^- (115) 

The same wave function in the ^^-representation reads: 

^loMk') = -^{-f'^{Sk)e-^-. (116) 

The wave function of the N = 2 state in the r- and ^-representations reads: 

^2o..(r) = {^Y^\-r^\Sf)il - |ar^)e-t^^ , (117) 

6 TT 5 

^2oMk') = ^(-)^/^^(5T)(-5 + 2^)e-ti . 
Vl5 a Ja a 



The mean radius squared of the 2 PoQ'Q'-state is equal to 

<'^'>23Po=^. (118) 

Mean radii squared, eqs. ( |115|) and ( |118| ), can be used for fixing the wave functions of the 
^0 and 23Po 



l^Pn and 23Pn states. 



6 D-Matrix Amplitude and K-Matrix Representation 

In this Section we consider the connection between the D-matrix amplitude and the 
amplitude in the K-matrix representation. Usually the analysis of experimental data is 
performed in the framework of the K-matrix amplitude: such an analysis has been carried 
out in refs. 0, ^. Performing a more complicated analysis, such as the D-matrix one, it 
is helpful to use the results of the K-matrix analysis at full scale. 

In Section 6.1 the transition from D-matrix amplitude to the K-matrix representation 
is considered for the case of a separate resonance decaying into a number channels. In 
Section 6.2 the case of two overlapping resonances is investigated. The results of the 
K-matrix analysis of the 00++ wave obtained in ref.p|, are quoted in Section 6.3. These 
results are used in the next Section as a guideline for the construction of the D-matrix 
amplitude. 
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6.1 One- Resonance Amplitude 

The multichannel amplitude for a separate resonance is given by eq.(0). The transition 
from this amplitude to the K-matrix representation is rather simple: 



\J Pn[s)Ann'[s)\l Pn'[s) - - 



n" n" 

(119) 
where 

g(n) (s)g(^') (s) JpJs)pn'(s) 

Let us ephasize that the D-matrix amplitude gives factorized K-matrix elements that 
means, for example, KnnKn'n' = Knn'Kn'n- This specific property results in a simplified 



K-matrix presentation of the amplitude what is seen at eq. ( |119|) . The K-matrix pole is 
defined by the equation: 

m^-s-ReB{s) = , 5(s) = ^ ^(""^(s). (121) 

n" 

Suppose that this equation is satisfied at some point s = p'^. Then in the vicinity of this 
K-matrix pole the K-matrix elements are of the form: 

This is exactly the form used in refs.p, §] for fitting experimental data. The pole position 
of the K-matrix elements is shifted, as compared to the propagator pole, in the magnitude 
ReB{p'^). Position of the K-matrix pole also differs significantly from the position of the 
amplitude pole, provided the value ImB{s) is not suppressed. This imaginary part should 
be saturated by meson-meson states, therefore their contribution can be considered as an 
additional "dressing" of bound state by a cloud of real mesons. Because of that we call 
K-matrix poles, which do not include this dressing, as poles of "bare" states P, p[. 



Another important consequence of eq. (|122|) is that the K-matrix coupling constants 



are connected with the D-matrix couplings as follows: 

g^''\K - matrix) = , ^^ ' = . (123) 

^1 + ReB'ip^) 

Therefore, if the D-matrix couplings satisfy quark combinatoric relations, the same re- 
lations are valid for the K-matrix coupling constants: the difference between g^"'\K — 
matrix) and g^"'\p'^) lays only in the factor (1 -|- ReB'{s))^^'^ which is common for all 
channels. 

6.2 Two-Resonance/Two Channel Amplitude 

Now we consider the case when two overlapping resonances decay into two meson-meson 
channels. Being rather simple for the investigation, this example contains all principal 
points of a representation of the D-matrix amplitude in the K-matrix form. 
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The D- matrix for the case under discussion is given by eq.(0) while the transition 
amphtude is written in eg. (pTf) . For the scattering amphtude in the channel 1 one has: 



{Ml - s)iMi - s) - BuB2i 



Pi{s)An{s) = — ^ ^^7^^72 ^^ ^-^ = (124) 



Ku+t{Ki2K2i-KnK: 
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1 - t{Kn + K22) + {K^2K2i - K^iK22) ' 

Remind that Mf = im?^ — Baa{s). The main question which should be answered here is: 
under what condition the quark combinatoric rules are valid for both D-matrix and K- 
matrix coupling constants? For this purpose, consider in more detail the matrix elements 
Kii and 7^22- The K-matrix form of eq.( |124D gives us the following expressions for Kn 
and K22'- 

^ 9f\^\ - s)^ + ^rVf - s)gt^ + 2ReB,29'?9't .,,.. 

{^,l{s) - s){^,l{s) - s) - {ReB,2? ' ^ ^ 

where 

lil{s) =ml- ReBUs) =ml-Y. ReB^^^Ks) . (126) 

n 

Below we show that the K-matrix couplings sustain the quark combinatoric rules if ReBi2 
is small . Indeed, in this case one has: 

^ ^i9!%)l^ii%)l .127) 

/if (s) - S /i^(s) - S 



An expansion of the right-hand side of eq. ( |127|) near the K-matrix poles gives: 



{qi [K — matrix)Y {§2 \K — matrix))" 



-'^nn — 2 - . 5/ _\ _ ' J'n-n i \^^^) 



where 



gi''\K - matrix) = , ^" ^^"^ =, (129) 



so the K-matrix couplings g^^' (K — matrix) satisfy the same relations as 9a if^a) because 



they differ by the factor Jl + ReB'^^dxD which is common for all channels. The K-matrix 
elements in the form given by eq. ([T0|) were used in refs.0, § for fitting experimental 
data. 

One sees that using quark combinatorics for the K-matrix couplings is justified if 
ReB 12 is small or, more generally, if the real parts of the non-diagonal loop transition 
diagrams are small. In the case of gg-states, a suppression of the non-diagonal terms 
ReBab{s) occurs because of orthogonality of {qq)a and {qq)b wave functions. For mesons 
— members of the same nonet — the suppression happens due to cancellation of the nn- 
and ss-loop diagrams (this cancellation is precise, provided the flavour SU(3)-symmetry 
is not broken). For mesons which are members of different nonets, the suppression of 
ReBab{s) is the result of orthogonality of radial wave functions. 
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The situation with imaginary parts of non-diagonal loop diagrams may be quite dif- 
ferent: the imaginary parts are determined by decays into channels with real mesons, and 
therefore they strongly depend on the location of thresholds. For example, two lowest K- 
matrix bare states, which are located below 1.2 GeV [Q, ||, despite of their being members 
of the same nonet, strongly mix due to large imaginary part of the loop transition dia- 
gram. The reason is that for the low-mass region only imaginary part of the non-strange 
quark loop diagram is properly saturated by tttt mesons, while the imaginary part of ss 
loop diagram, which should be saturated by KK, rjr] and rjrj' mesons, is small because of 
the threshold suppression of phase spaces, and it is not able to cancel the tttt contribution. 

Another situation with the D-matrix and K-matrix couplings may happen when qq 
states mix with a particle of another origin, for instance, with a glueball. This case is a 
subject of our further investigation. 

6.3 Results of K-Matrix Fit of 00++ Amplitude 

Analysis of scalar/isoscalar states in the mass region below 1900 MeV was performed in the 
framework of the K-matrix approach in ref. basing on the following data set: GAMS 
data for n^p — > TT^ir^n |2^, rir]n ||2l|, rjrj'n ||2^; CERN-Miinich data for Tc^p -^ n^Ti^n 



[^ ; Crystal Barrel data for pp -^ tt^tt^tt'^, vt^tt^?], ti^t]-)] p4|, |2^; BNL data for tttt — > KK 
[ p^ . The K-matrix elements were used in the form similar to eq. (127): 

7^ sr^g^aKK -matrix) g^^'\K -matrix) 

J^nn' = Z^ 5-3 \- Inn' , (loOj 

with five channels taken into account n = tttt, KK, r]r], r]r]', rnnm. Couplings g^\K — 
matrix) are supposed to obey quark combinatoric relations that allows to reconstruct the 
gg/glueball content of bare states below 1900 MeV. The simultaneous fit of the two-meson 
spectra of refs. [pO| j-p6| gives two solutions which describe well the data sets: 



Solution I. Two bare states f^''''''{720 ± 100) and /o^'^'^^(1260 ± 30) are members of the 
l^Po qq-nonet, with /q°^'^(720) being ss-rich state, 0(720) = — 69°± 12° (the flavour wave 
function is defined as nnCos(j)+ssSin(j) ). The bare states /g"'"^(1600±50) and /q°^*^(1810± 
30) are members of the 2^Po-nonet; /o'"'^(1600) is dominantly nn-state: 0(1600) = — 6° ± 
15°. The state /q°^^(1235 ± 50) is superfluous for gg-classification, being a candidate 
for the lightest glueball: its couplings to two-meson channels obey quark combinatoric 
relations for glueball. 
Concluding, we have the following states in the solution I: 

Type of state Flavour wave function 

/o^»'''^(720) -^ l^Poqq 0.37nn - 0.93ss 

/o^"''^(1260) -^ l^Poqq 0.93nn + 0.37ss , . 

/o^"^'"(1600) ^ 23Pogg 0.995nn - O.lOss ^^ 

/o^"^'"(1810) -^ 2^Poqq O.lOnn + 0.995ss 

f^"'''^{1235) -^ Glueball -^ OMnn + 0A2ss . 

For the glueball, "flavour wave function" refers the qq intermediate state in the decay 
process (see fig. 6c). 
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Solution II. Basic nonet members are tlie same as in solution 1. The members of 
the 2^Po-nonet are the following: /q"^'^(1235) and /o"'''^(1810); both these states have 
significant ss-component: 0(1235) = 42° ± 10" and 0(1810) = -53° ± 10°. The state 
fbare^-^^QQ _|_ gg-j jg superfluous for the gg-classification, being a good candidate for the 
lightest glueball. 

So, we have in the solution II the following set of states: 

Type of state Flavour wave function 

/o*»'''^(720) -^ l^Poqq 0.37nn - OMss 

/^""-^UeO) -^ I'^Poqq 0.93nn + 0.37ss , . 

ybarej^235) ^ 2^Poqq 0.74nn + 0.67ss ^ ' 

/o^"^'^(1810) -^ 2^Poqq 0.67nn - 0.74ss 

/o*"''"(1560) -^ Glueball -^ 0.91nn + 0A2ss . 

For both solutions five scalar resonances with nearly coinciding positions of the amplitude 
poles were found: 

Resonance Position of pole on the 

complex-M plane, in MeV units 

/o(980) 1015±15-i(43±8) 

/o(1300) 1300±20-i(120±20) (133) 

/o(1500) 1499±8-z(65±10) 

/o(1750) 1750±30-i(125±70) 

/o(1200/1600) 15301^^0 - ^(560 ± 140) . 

Comparison of the positions of bare states and those of resonances shows that the observed 
states are strong mixtures of the bare states and two/four-meson states into which these 
resonances decay. 

Concluding this section, let us discuss the problem of applying quark combinatorics 
to the K-matrix decay couplings for the bare states in the mass region 1200-1600 MeV, 
/o^'^'"^(1260), /o*'^^'=(1235) and /o*'^'''=(1600): one of them is supposed to be the glueball. This 
means that corresponding non-diagonal loop diagrams for the transitions Glueball -^ 
l^Poqq and Glueball -^ 2^PQqq are not small, and the arguments of Section 6.2 do not 
work. Nevertheless, eqs. (131) and (132) show us that quark combinaroric rules can 
be applied to these resonances as well. The reason is that the discussed gg-states have 
approximately the same nn/ss content as the intermediate state in the glueball decay 
(the flavour wave function of this intermediate state is written in eqs. (131) and (132)). 
In this case, it is obvious that quark combinatoric relations are almost the same both for 
the K- and D-matrix amplitudes, see eq. ( |125|) . 

7 D-Matrix Fit of the 00++ Amplitude 

In this Section we consider the mixing of the glueball state with neighbouring gg-states. 
We restrict our consideration by the region 1100-1900 MeV: the calculation of the lowest 
state /o(980) requires to take into account the mechanism which saturates the imaginary 
part of the quark-antiquark loop diagram by meson-meson states in the precise form. In 
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the large-mass region, above the main meson-meson thresholds, the detailed description 
of the saturation is not important. 

7.1 Parameters and Results 

For the calculation of the quark loop diagrams, we should fix the vertices ga{s). We 
parametrize the vertices for the transition state a -^ nn in the simple form: 

l^Po qq - state : gi{s) = 71 ./i |i^; 

2^Po qq — first state : (72 (s) = 72 \^ 



d-. 



fc2+o-2 k'^+cr2+h 



(134) 



Gluehall : g,{s) = 73 ^ fg; 

2'^Po qq — second state : gi{s) = g2{s). 

2 

Here fc^ = | — rri^ and k"^ = ^ — rn^ where m is the constituent quark mass (fixed at 
350 MeV for non-strange quark and at 500 MeV for strange one); rria, 7a and Ca are 
parameters. Factor d is due to orthogonality of the l^Poqq and 2^Poqq states: we put 
Re Bi2{so) = at ^/sq = 1.5 GeV. (In the case of the s-dependent P-functions the 
orthogonality requirement for loop transition diagrams cannot be fixed at all values of s) . 

The parameters m^, o"a, h and 7^ (a = 1, 2, 3) should be determined by masses and 
widths of the physical resonances. However, the masses rua are approximately fixed 
by the K-matrix fit of ref. 0, where masses of the K-matrix poles, /i^'^'''^, are defined: 
/i^ ~ ml — Bl^{jj?g) (see also eqs. (131) and (132)). Let us stress that m^ is the mass of 
pure gluonic glueball which is a subject of Lattice QCD calculation. 

Parameters which are found in our fit of the OO"*"^ amplitude in the mass region 1200- 
1900 MeV are given in Table 3 for the solutions I and II. Using these parameters, we 
calculate the couplings aa which are introduced by eqs. (^7|) and (|38|): 



|/o(1300) >^ |/o(1500) >^ |/o(1530) >^ |/o(1780) >^ . (135) 

These couplings determine relative weight of the initial state a in the physical resonance 
A: 

Wa{A) = la.p (136) 

The probabilities Wa are given in Table 3 together with masses of physical resonances. 
Ma, and masses of input states, nia- 

7.2 Glueball/gg-State Mixing 

In order to analyze the dynamics of the glueball/gg mixing, we use the following method: 
in final formulae the vertices are replaced in a way: 

ga{s)-^iga{s), (137) 
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with a factor ^ running in the interval < ^ < 1. The case ^ = corresponds to switching 
off the mixing of input states. Input states are stable in this case, and corresponding poles 
of the amplitude are at Sa = Tn\. Fig. 8 shows the pole position at .^ = for solution I 
(fig. 8a) and solution II (fig. 8b). For the glueball state, m^ is the mass of a pure glueball, 
without qq degrees of freedom. In solution I the pure-glueball mass is equal to 

rripure giuebaii{Solution I) = 1225 MeV, (138) 

that definitely disagrees with Lattice-Gluodynamics calculations for the lightest glueball. 
In solution II 

rripure giuebaiiiSolution H) = 1633 MeV. (139) 

This value is in a good agreement with recent Lattice-Gluodynamics results: 
1570±85(stat)±100(syst) MeV |, |] and 1707±64 MeV §. With increasing ^ the poles 
are shifted into lower part of the complex mass plane. Let us discuss in detail the solution 
II which is consistent with Lattice result. 

At ^ ~ 0.1 — 0.5 the glueball state of solution II is mixing mainly with 2^Po qq-state, 
at ,^ ~ 0.8 — 1.0 the mixture with I^Pq qq-sta.te becomes significant. As a result, the 
descendant of the pure glueball state has the mass M = 1450 — z450 MeV. Its gluonic 
content is 47% (see Table 3). We should emphasize: the definition of Wa suggests that 
J2A=i,2,3,4^giuebaii{,A) ^ 1 because of the s-dependent Bah in the propagator matrix. 

Hypothesis that the lightest scalar glueball is strongly mixed with neighbouring qq 
states was discussed previously (see refs. 0, p^; &nd references therein). However, the 



attempts to reproduce a quantitative picture of the glueball/gg-state mixing and the mass 
shifts caused by this mixing cannot be successful within standard quantum mechanics ap- 
proach that misses two phenomena: 

(i) Glueball/gg-state mixing described by the propagator matrix can give both a repulsion 
of the mixed levels, as in the standard quantum mechanics, and an attraction of them. 
The latter effect may happen because the loop diagrams B^h are complex magnitudes, 
and the imaginary parts ImBab are rather large in the region 1500 MeV. 
(ii) Overlapping resonances yield a repulsion of the amplitude pole positions along imaginary- 
s axis. In the case of full overlap of two resonances the width of one state tends to zero, 
while the width of the second state tends to the sum of the widths of initial states, 
T first — and Ts^cond — Fi + F2. For three overlapping resonances the widths of two 
states tend to zero, while the width of the third state accumulates the widths of all initial 
resonances, Tmrd ^ Fi + F2 + F3. 

Therefore, in the case of nearly overlapping resonances, what occurs in the region near 
1500 MeV, it is inevitable to have one resonance with a large width. It is also natural that 
it is the glueball descendant with large width: the glueball mixes with the neighbouring 
l^Po QQ and 2^Po qq states, which are both nn rich, without suppression. 

8 Conclusion 

On the basis of N/D method we have developed the matrix propagator approach (D- 
matrix method) for the analysis of hadron spectra in the case of overlapping resonances. 
This technique preserves the unitarity of a scattering amphtude as well as its analyticity 
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in the whole complex-s plane. Being rather similar to the K-matrix technique for the 
case of overlapping qq resonances the D-matrix technique can give different results when 
a particle of another nature (e.g. glueball or hybrid) interacts with a qq resonance. The 
case of the interaction of the scalar glueball with neighbouring qq states in the region 1100- 
1800 MeV is considered on the language of quarks and gluons. The content of the initial 
qq and glueball states in the physical resonances was calculated. It appeared that the 
lightest gluodynamic glueball is dispersed over neighbouring resonances mixing mainly 
to I^Pq qq and 2^Po QQ states. With this mixing the glueball descendant transforms 
into broad resonance, /o(1530l25o)- This resonance contains (40-50)% of the glueball 
component. Another part of the glueball component is shared between comparatively 
narrow resonances, /o(1300) and /o(1500) which are descendants of l^Po QQ and 2^Pq qq 
states. 

The important development for understanding of meson spectra would be the direct 
application of the D-matrix approach (instead of the K-matrix one) to the analysis of 
the meson production amplitude. Such approach will give not only the correct position 
and mixing of the exotic states with the gg-mesons but also clarify the situation at low 
energies where a usual K-matrix approach starts to violate analytical properties of the 
scattering amphtude. 
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Table 1 

Coupling constants given by quark combinatorics for a glueball decaying into two 

pseudoscalar mesons in the leading and next-to- leading terms of the 1/N expansion. $ 

is the mixing angle for nn and ss states, and O is the mixing angle for rj — rj' mesons: 

1] = nn cos Q — ss sin O and rj' = nn sin + ss cos 0. 
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Table 2 

Coupling constants given by quark combinatorics for a (gg)a-nieson decaying into two 

pseudoscalar mesons; here {qq)a = nn cos ^ + ss sin $. 



Channel 


The (gg)a-meson decay 

couplings in the 
leading terms of l/N^ 


The (gg)a-meson decay 
couplings in the next- 
to-leading terms of 1/Nc 
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expansion (Fig. 7b) 
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Table 3 
Masses of the initial states, coupling constants and qq/ ghiehaW content of physical states. 
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Fig. 1. c) Meson states in the infinitely increasing quark model potential b) 
Conventional picture for the transition qq — states -^ meson states. 
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B„>^ +^=4 B„>==^ B, 






rrii XX rn, x x m 





d) 



Fig. 2. Diagrams responsible for the decay /mixing process. 
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e) 



Diagrams for entirely-composite particle. Fig. 3. 
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Fig. 4. Stable (a) and non-stable (b) levels in the potential model. 
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q) 






q 

c) 



e) 



Fig. 5. The gg-interaction diagrams at r < Rconfinement (a) presented as a set of pole 
terms (b). Saturation of the gg-loop diagram at r > Rconfinement (c) by a set of meson 

loop diagrams (d,e). 
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Fig. 6. Diagrams responsible for glueball decay. 
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meson 




Fig. 7. Diagrams of (gg)^- meson decay. 
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Fig. 8. Complex-y^ plane (M = Re^/s, — r/2 = Im^/s): location of 00++ amplitude 

poles after replacing ga -^ ^Qa- The case ^ = gives the positions of masses of input qq 

states and gluodynamic glueball; ^ = 1 corresponds to the real case. 
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